Properties of asymmetric nuclear matter in different approaches 
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Properties of asymmetric nuclear matter are derived from various many-body approaches. This 
includes phenomenological ones like the Skyrme Hartree-Fock and relativistic mean field approaches, 
which are adjusted to fit properties of nuclei, as well as more microscopic attempts like the Brueckner- 
Hartree-Fock approximation, a self-consistent Greens function method and the so-called Vi om u ap- 
proach, which are based on realistic nucleon-nucleon interactions which reproduce the nucleon- 
nucleon phase shifts. These microscopic approaches are supplemented by a density-dependent con- 
tact interaction to achieve the empirical saturation property of symmetric nuclear matter. The pre- 
dictions of all these approaches are discussed for nuclear matter at high densities in /3-equilibrium. 
Special attention is paid to behavior of the isovector component of the effective mass in neutron-rich 
matter. 
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I. INTRODUCTION 



The equation of state for nuclear matter under exotic conditions is one of the main topics in modern nuclear physics. 
This interest in the properties of nuclear matter at high densities and large proton-neutron asymmetries is partly 
motivated by the fact that this information is required in theoretical models of compact objects in astrophysics like 
neutron stars or the simulation of supernovae. However, the study of nuclear systems with large isospin asymmetries 
are also the subject the forthcoming radioactive-ion-beam facilities such as FAIR at GSI or SPIRAL2 at GANIL. 

Different ways have been developed to obtain predictions for the properties of nuclear systems under exotic condi- 
tions. One way is to start from phenomenological models which successfully describe the properties of stable nuclei. A 
very popular approach along this line is the use of an effective density dependent Skyrme- type interaction[l|,[2,[3,|3,[l]. 
Modern Skyrme parameterizations have been developed, which were constrained in their fitting procedures to obtain 
results for neutron-rich nuclear matter compatible to those of microscopic calculations. Here we mention the Skyrme- 



' Lyon (SLy) forces @, which have been used in studies of the neutron-star crust @, Q. 
£T) . Also the relativistic mean-field approximation has very successfully been used to describe the properties of stable 
nuclei . Attempts have been made to derive a set of density-dependent meson- nucleon coupling constants of density 
I dependent relativistic mean field (DDRMF) calculations from microscopic Dirac Brueckner Hartree Fock calculations 
with a readjustment in order to reproduce the bulk properties of stable nuclei and the saturation point of symmetric 
OO | nuclear matter [Io|. 

The parameters of such Skyrme Hartree Fock and DDRMF calculations have been fitted to the data of stable nuclei 
and the predictive power of these simple phenomenological approaches could be rather limited. This would mean that 
the predictions for nuclear systems with exotic values for density and proton-neutron asymmetries may not be very 
reliable. 

Therefore we will also consider some so-called microscopic approaches, which start from models of the nucleon- 
nucleon (NN) interaction, which are adjusted to describe the experimental phase shifts of NN scattering at energies 
below the pion threshold [TJ The traditional models of such realistic NN interactions like e.g. the charge- 

dependent Bonn potential CDBONNfll] or the Argonne potential V18[l2j contain rather strong short range compo- 
nents, which make it inevitable to employ non-perturbative approximation schemes for the solution of the many-body 
problem for the nuclear hamiltonian based on such interactions [l3j . Such non-perturbative approximations include 
the Brueckner hole-line expansion with the Brueckner Hartree Fock (BHF) approximation, the self-consistent evalua- 
tion of Green's function using the T-matrix approximation [tH fl5L ITo, []~7L Il8l. Ixol. I20I |2~0] (SCGF) and also variational 
approaches using correlated basis functions [23, [H, [24| . 

If such microscopic calculations would reproduce the properties of nuclear systems under normal conditions, one 
could argue that a scheme which reproduce the data of two nucleons in the vacuum (NN phase shifts) as well as nuclear 
data at normal densities should also provide reliable results for nuclear systems at densities beyond the saturation 
density of nuclear matter. Unfortunately, however, such microscopic calculations fail to reproduce the saturation point 
of symmetric nuclear matter or the bulk properties of finite nuclei with good precision [r|. For example in nuclear 
matter such calculations yield results for the saturation point, which are located on the so-called Coester band(25[, 
i.e. they either yield too little binding energy or a saturation density well above the empirical value of po = 0.16 fm -3 . 
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In recent years large progress has been made developing tools for essentially exact calculations for nuclear system 
with mass number up to about A = 12 (26|. These calculations demonstrate that using realistic models for the NN 
interaction in a non-relativistic hamiltonian precise results for nuclear systems are obtained if and only if these realistic 
two-body interactions are supplemented by a three-body force. 

Therefore we will consider in the present study the results of BHF and SCGF calculations, adding a simple density- 
dependent contact interaction which is adjusted to describe the saturation point of symmetric nuclear matter. 

During the last years a different calculation scheme has evolved, which is also based on the NN scattering data but 
tries to decouple the low- and high- momentum components in the nuclear hamiltonian using renormalization group 
methods]!?], [H, [29|, [3(j. The interaction resulting for the low-momentum regime of nuclear structure, which we will 
refer to as V[ ow k is rather soft, which implies that non-perturbative tools of many-body theory can be applied_3lj. Also 
it is very attractive, that Vi ow h turns out to be essentially model-independent, if the cut-off for the low-momentum 
regime is chosen appropriately. 

However, evaluating the properties of nuclear matter by using Vi ow k in a Hartree-Fock or BHF calculation, one 
does not obtain a saturation point [HI, the binding energy per nucleon increases with increasing density and 
nuclear systems tend to collapse to high densities. This will be compensated by adding a density-dependent contact 
interaction, which is adjusted in the same fashion as for the BHF and SCGF cases discussed above. 

Therefore we will apply two different schemes: two phenomenological methods (Skyrme Hartree Fock and relativistic 
mean field DDRMF) employing parameterizations, which tend to reproduce some features of microscopic calculations 
and three microscopic approaches (BHF, SCGF and Hartree-Fock with Vi ow h)i which are supplemented by a simple 
contact interaction to reproduce the empirical saturation properties of symmetric nuclear matter. For all these methods 
we study bulk features of nuclear matter at large densities and proton neutron asymmetries. A comparison of the 
resulting values for nuclear compressibility, symmetry energy, proton fraction in (3 equilibrium and effective masses, 
which characterize the density of states around the Fermi energy, should provide information about the uncertainties 
in the extrapolation of nuclear properties in exotic regions of density and asymmetry. Also we expect some hints 
about the reliability and problems of the individual approaches. 

After this introduction the section two shall briefly outline the different approximation schemes which we are going 
to consider. The results are discussed in section three and the main conclusions are summarized in the final section. 



II. DIFFERENT MANY-BODY APPROACHES 



A. Skyrme Hartree-Fock 



A very popular many-body approach in nuclear physics is the Skyrme-Hartree-Fock approach which can be found 
e.g. in [3, 0, H, 13] • The Skyrme interaction leads to an energy functional 



E = J H(r)dr, (1) 



where H is the Hamiltonian densityin the Hartree-Fock approximation. In case of infinite asymmetric nuclear matter 
the Hamiltonian density writes [2j, [6| 

H = Hk + H c ft + Ho +H 3 (2) 

where Hk is the kinetic energy term, H e s an effective mass term, Ho a zero range term and H 3 a density dependent 
term. These terms are given by 

H cS = ±[il(2 + .T 1 )+i 2 (2 + .T 2 )]T ( 

+§ [t 2 (2x 2 + 1) - t 1 {2x 1 + 1)] [t pPp + r nPn ] , 
Ho = \t [(2 + x )p 2 ~(2x + l)(p 2 p + pl)], 

n 3 = ±t 3 p a [(2 + x 3 )p 2 -(2x 3 + l)(p 2 p + pl)}, (3) 

where the coefficients ti, Xi, Wo, and a are the parameters of a generalized Skyrme force. In the present study we 
apply the commonly used parameterization SLy4. 
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The densities p and t are defined in terms of the corresponding densities for protons and neutrons p — p p + p n and 
t = T p + r n . If we identify the isospin label (i = p,n), the matter densities for protons and neutrons are given by 

„ f kFi d 3 k 1 - 
P > = 2 J (2^ = 3^' (4) 

where fcpi denotes the corresponding Fermi momentum and the spin degeneracy gives a factor of 2. We obtain the 
kinetic energy density by 



n = 2/ 7^ 2 = A(3^) 5/3 . (5) 



d 3 k 2 _ 1 

(2^)3 ~ 57 2 

The Hamiltonian density leads to the energy per nucleon E/A = TL/ p and the single particle energy may be written 
as a function of momentum k 

«<*)- ^ + v, (0) 

with an effective mass m* and a density dependent Potential VI. 



B. Brueckner-Hartree-Fock 



Starting from realistic NN interactions, we have to use more advanced many-body approximations like the 
Brueckner-Hartree-Fock (BHF) approach, which have the capability to account for the effects of correlations, which 
are due to the strong tensor and short-range components of such realistic NN interactions. In the BHF approximation 
this is achieved by evaluating the so-called G-matrix, which corresponds to the in-medium scattering matrix. The 
self-energy of a nucleon with isospin i, momentum k and energy to in asymmetric nuclear matter is defined in the 
BHF approximation by [HI, H(| 

^f HF = E <*9|G(n)|kg) y nS(g). (7) 

3 

In this equation n°(g) refers to the occupation probability of a free Fermi gas of protons (j — p) and neutrons 
(j = n) like in the mean-field or Hartree-Fock approach. This means that for asymmetric matter with a total density 
p = p p + p n this probability is defined by 

n > {q) -\0 io r \q\>k F3 , (8) 

with Fermi momenta for protons (kp p ) and neutrons (kp n ). 

The antisymmetrized G matrix elements in eq. ([7]) are obtained from a given NN interaction by solving the Bethe- 
Goldstone equation 

(fcg|G(fi)|fcg) tf =(kq\V\kq) i:j + J d 3 Pl d 3 p 2 (kq\V\ Pl p 2 ) l: j 

x Q(PlhP2j) (9) 

Q - (e Pu i + s P2tj ) + irj 
x (pip 2 |G(0)|fcg)y. 

The single-particle energies e P i of the intermediate states should be the corresponding BHF single-particle energies 
which are defined in terms of the real part of the BHF self-energy of eq. ([7]) by 

E kl = ^+Rc[^f HF {k,uj = s kl )\, (10) 
2m 

with a starting energy parameter £1 = ui + e q j in the Bcthc-Goldstone equation (J9j> . 
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The Pauli operator Q(pii,P2j) restricts the intermediate states to particle states with momenta p\, p2, which are 
above the corresponding Fermi momentum. However, the single-particle spectrum is often parameterized in the form 
of an effective mass 

k 2 

e*i » TT^T + ^ ( n ) 
2m* 

so that a so-called angle-averaged propagator can be defined, which reduces the Bethe-Goldstone equation to an 
integral equation in one dimension. The exact Pauli operator has been treated in [33| ■ 



C. Relativistic models 



An interesting extension of the BHF approach is the Dirac-Brueckner-Hartree-Fock (DBHF) approach which ac- 
counts for relativistic effects as well as correlation effects as described in terms of the G matrix. In the DBHF 
approximation one evaluates the self-energy Sj(fc) in a way very similar to the BHF approximation keeping track of 
the Lorentz structure. Attempts have been made to fit a density dependent Relativistic Mean Field model (DDRMF) 
to the results of such DBHF calculations [l(| [H[ • Here we employ such a DDRMF model which has been obtained by 
fitting density dependent coupling constants for the meson-nucleon vertices to reproduce the self-energy of the DBHF 
calculations of van Dalen et al. [35J, which was based on the Bonn A potential. 

Both approaches, the DDRMF as well as the DBHF, start from a Lagrangian, which includes baryons, mesons and 
the interaction 

C = C B +C M + C int . (12) 



However, they differ in the included mesons and the coupling operators [UJ, [35|, [36|, l37( . In the DDRMF model the a, 
5, uj, and p mesons are included in scalar and vector couplings, respectively. The variation of the Lagrangian leads to 
a Dirac equation for the nucleons 

[7 • k* + m* Di ] u(k, s, i) = 70 E* u(k, s, i), (13) 
where u(k, s, i) denotes the plane wave Dirac spinor with momentum k, spin s and isospin i 




u(k,s,i) = J 1 2 ^ Dt \ ) xi/ 2 (s)xi/ 2 (j). ( Ul) 

The starred quantities contain the different components of the nucleon self-energy: the scalar, time-like vector and 
space-like contributions 

m* Di = M + Z Sji (k,k F ) 
E* = Ej(fc) - £„,*(*,**■) 

k* = ki + kiT, Vt i(k,k F ), (15) 

where fej is the unit vector along the momentum ki of the nucleon. The general form of the nucleon self-energy in 
infinite matter is obtained by evaluation of the meson exchange in spin saturated nuclear matter 

Si(fc, k F ) = S s ,i(fc, k F ) + 70 £o,i(fc, k F ) +7 • ki £v>(&, k F ). (16) 

In the fit ting process special attention has to be payed to the rearrangement contribution to the time- like vector 
self-energy [381 ] . the relativistic operator structure [391 ]. and the proper renormalization due to the space-like vector 
contribution to the self-energy |l0| . 

The single particle energy in the DDRMF model is obtained from the Dirac equation |T 



Ei{k) = Vfef +m^ 2 + £ 0li! (17) 



and the energy-momentum tensor leads to the energy density in asymmetric nuclear matter 

" Jo 

i—p.n w 



2 

i— p,n 
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with the scalar densities p\ and the baryon densities pi 



k 2 dk. (18) 



D. Self-consistent Green's function 



One of the drawbacks of the BHF approximation is the fact that it does not provide results for the equation of state, 
which are consistent from the point of view of thermodynamics. As an example we mention that BHF results do not 
fulfill e.g. the Hugenholtz van Hove theorem. This is due to the fact that the BHF approximation does not consider 
the propagation of particle and hole states on equal footing. An extension of the BHF approximation, which obeys 
this symmetry is the self-consistent Green's function (SCGF) method using the so-called T-matrix approximation. 
During the last years techniques have been developed, which allow to evaluate the solution of the SCGF equations 
for microscopic NN interactions [lH 

HE E3, EH, Those calculation demonstrate that for the case of realistic 
NN interactions, the contribution of particle-particle ladders dominates the contribution of corresponding hole- hole 
propagation terms. This justifies the use of the BHF app roximation and a procedure, which goes beyond BHF and 
accounts for hole-hole terms in a perturbative wav[iol l4ll]. This leads to a modification of the self-energy in the BHF 
approximation by adding a hole-hole term of the form 



A£f lp (fc,c)=]T / d 3 p / d 3 hj 

„• Jk Fj JO JO 



FS d 3 h x J d 3 h 2 

(19) 



j 
x 

lu + E P j - e hl i - £h 2 j - if] 
The quasi-particle energy for the extended self-energy can be defined as 

€i = £ + MZ? HF (k,u = e%) + AE? w *(fc )W = eZ)l (20) 

Accordingly, the Fermi energy is obtained evaluating this definition at the Fermi momentum k = kpi for protons and 
neutrons, respectively, 

£Fi = eZ i . (21) 

The spectral functions for hole and particle strength, SMk, ui) and Sf(k, oS), are obtained from the real and imaginary 
part of the self-energy S = T, BHF + AS 2 "? 

S h(p) (k uj) = ±- ImSjfoq;) 

1 K ' 1 7r[w-fc 2 /2m-Re£,(fc,w)] 2 + [Im£ 4 (fc,w)] 2 ' V 1 

where the plus and minus sign on the left-hand side of this equation refers to the case of hole (h, ui < epi) and 
particle states (p, lo > e_Fi), respectively. The hole strength represents the probability that a nucleon with isospin i, 
momentum k, and energy u> can be removed from the ground state of the nuclear system with the removal energy u>, 
whereas the particle strength denotes the probability that such a nucleon can be added to the ground state of the 
system with A nucleons resulting in a state of the A + 1 particle system which has an energy of lu relative to the 
ground state of the A particle system. Hence the occupation probability is obtained by integrating the hole part of 
the spectral function 

m(k) = f F * duS^(k,u). (23) 



Note that this yields values for the occupation probability, which ranges between values of and 1 for all momenta 
k, leading to a partial depletion of the hole-states in the Fermi gas model (k, < kp) and partial occupations for states 
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with momenta h > kp. A similar integral yields the mean energy for the distribution of the hole and particle strength, 
respectively 

f eFz duuSHk,u) 

< £ **» - — „, w ■ < 24 > 

| du> u> Sf(k,u>) 

M*)> = ■ ( 25 ) 

Our self-consistent Green's function calculation is defined by identifying the single particle energy in the Bethe- 
Goldstone equation as well as in the 2/ilp correction term in eq. (|19|) by 



| (Shrik)) for k < kp T , 
kT \( £ P T(k)) ioxk>k Fr . 

This definition leads to a single particle Greens function, which is defined for each momentum k by just one pole 
at to = £fc r . Hence in the calculation of the self-energy the mean value of the spectral distribution is considered. 
However, the modified occupation of nucleons obtained by the spectral functions are not included in the calculation 
of the self-energies. The total energy per nucleon is evaluated by 

E gj I d 3 k JJ2 du SHk, u;)(k 2 /2m + u)/2 

A E t Jd 3 kn t (k) ■ ( > 



E. Renormalization of the NN interaction 



It is very reasonable to assume that the long-range or low-momentum part of the NN interaction is fairly well 
described in terms of meson-exchange, while different (quark) degrees of freedom arc getting important to describe 
the short-rang or high-momentum components of the NN interaction. Therefore it is quite attractive to disentangle 
these low-momentum and high-momentum components from each other. This means that one tries to define a model 
space, which accounts for the low-momentum degrees of freedom and renormalizes the effective hamiltonian for this 
low-momentum regime to account for the effects of the high-momentum parts, which are integrated out. 

This concept of a model space and effective operators appropriately renormalized for this model space has a long 
history in approaches to the nuclear many-body physics. As an example we mention the effort to evaluate effective 
operators to be used in Hamiltonian diagonalization calculations of finite nuclei. For a review on this topic see e.g. [42j . 
The concept of a model space for the study of infinite nuclear matter was used e.g. by Kuo et al. [iH H. |45| . 

During the last years this concept has received a lot of attention and led to the definition of the so-called Vi ow k 
interaction. One way to determine this interaction is to follow the unitary-model-operator approach (UMOA) [46] . 
We follow the usual notation and define a projection operator P, which projects onto the model space of two-nucleon 
wave functions with relative momenta k smaller than a chosen cut-off A. The operator projecting on the complement 
of this subspace is identified by Q and these operators satisfy the usual relations like P + Q = 1, P 2 = P, Q 2 = Q, and 
PQ = = QP. It is the aim of the Unitary Model Operator Approach (UMOA) to define a unitary transformation 
U in such a way, that the transformed Hamiltonian does not couple the P and Q space, i.e. 

QU- 1 HUP = 0. (28) 

The technique to determine this unitary transformation is nicely been described by Fujii et al.[47j (see also [3l|). It 
leads to an effective hamiltonian H e f f = ho + V e f / , which contains the term of the kinetic energy ho and an effective 
interaction V e f / given by 

V eff = V lowk = U- 1 (h + V)U-h . (29) 

Diagonalising this effective hamiltonian in the low-momentum model-space, one obtains eigenvalues which are identical 
to the diagonalisation of the original hamiltonian ho + V in the complete space. This means solving the Lipmann 
Schwinger equation for NN scattering using this Vi ow k with an cut-off A yields the same phase shifts as obtained 
for the realistic interaction V without a cutoff. One finds that the resulting Vi ow k 1S essentially independent of the 
underlying realistic interaction V, if is fitted to the experimental phase shifts and if the cut-off A is chosen around 
A = 2 fm , which means that the model space includes scattering up to the pion threshold. In that sense Vi ow k is 
unique and, as it reproduces the NN scattering phase shifts it can also be regarded as a realistic interaction like e.g. 
the CDBONN or Argonne V18 interactions. 
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FIG. 1: (Color online) Comparison of binding energy per nucleon of symmetric nuclear matter as obtained from 
Skyrme SLy4, DDRMF, BHF, SCGF, and Vlowk. Results of approaches based on realistic NN interactions are also 
compared with an additional contact interaction of the form displayed in ea l30l 

Since, however, the high-momentum or short-range components have been integrated out by means of the unitary 
transformation of eq. (j2"5)) , the Vi ow k does not induce any short-range correlations into the nuclear wave function. This 
leads to the nice feature that mean-field calculations using Vi ow k lead to reasonable results and corrections of many- 
body theories beyond mean field are weak[3lj]. On the other hand, however, it is this lack of short-range correlation 
effects, which are modified in the medium, which prevents the emergence of a saturation point in calculations of 
symmetric nuclear matter [28ll32j| . 

In order to to achieve saturation in nuclear matter one has to add three-body interaction terms or a density- 
dependent two-nucleon interaction. This is not very astonishing as it is known that a renormalization of a two-body 
operator leads to many-body terms [isl l49j|. Therefore it is quite natural to supplement the effective interaction Viowk 
by a simple contact interaction, which we have chosen following the notation of the Skyrme interaction to be of the 
form 

AH = \t QP 2 + ^Up 2+a , (30) 

where p is the matter density, to, t% and a are parameters. For a fixed value of a (typically a =0.5) we have fitted to 
and £3 in such a way that a Hartree-Fock calculation using Vi ow k plus the contact term of eq. (|30|) yields the empirical 
saturation point for symmetric nuclear matter. 

The same parameterization of a contact term has been used to evaluate corrections to the self-energy of BHF and 
SCGF in such a way that also these calculations reproduce the saturation of symmetric nuclear matter. 

Note that this contact term is an isoscalar term and does not influence the symmetry energy, proton fractions in 
/3-equilibrium, and effective masses. 

III. RESULTS AND DISCUSSION 

All results of calculations, which refer to realistic NN interactions, have been obtained using the CDBONN fi~l| 
interaction. This includes all BHF and SCGF calculations. Also the evaluation of Vi ow k has been based on the 
proton-neutron part of CDBONN. Using a cut-off parameter A = 2 fm -1 , these results do not significantly depend on 
the underlying interaction. The Skyrme Hartree-Fock calculations have been done using the parameterization SLy4 
and for the relativistic mean-field calculation the parameterization for DDRMF in [Uj has been used. 

First let us turn to the binding energy of symmetric nuclear matter, which are displayed in Fig. [T] Compared to 
other realistic NN interactions the CDBONN potential, which we have chosen here is a rather soft NN interaction 
with a weak tensor force. This is indicated by the results for the saturation point of symmetric nuclear matter as 
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SLy4 DDRMF BHF BHF(+ct) SCGF SCGF(+ct) V iowk + ct 


Po [fin" 3 ] 
E/A( Po ) [MeV] 
K [MeV] 
a s (po) [MeV] 


0.160 0.178 0.374 0.161 0.212 0.160 0.160 
-15.97 -16.25 -23.97 -16.01 -11.47 -16.06 -16.0 
230 337 286 214 203 270 258 
32.0 32.1 51.4 31.9 34.0 28.3 21.7 



TABLE I: Properties of symmetric nuclear matter are compared for Skyrme SLy4, DDRMF, BHF, SCGF, and Vi om k- The 
results, which are listed in the columns labeled with +ct are obtained employing the additional contact interaction of eq. (|30[) 
with parameters as listed in table HU The quantities listed include the saturation density po, the binding energy at saturation 
E/A, the compressibility modulus K and the symmetry energy at saturation density as(po)- 





BHF SCGF V lowk 


t [MeV fm 3 ] 
t 3 [MeV fm 3+3a ] 


-153 -311 -438.1 
2720 3670 6248 



TABLE II: Parameters to and t$ defining the contact interaction of eq. (|30p as obtained for the fit to the saturation point 
p = 0.16/m -3 and E/A = -16.0 MeV at a = 0.5 for various realistic approaches. 



obtained in the BHF approximation (the minimum of the dashed black line in Fig. [T] and data in table UJ ) . The 
saturation density is larger than twice the empirical value and the calculated energy is well below, which means that 
the CDBONN result is located in the large binding energy high density part of the Coester bandfl3l]. 

The hole-hole contributions to the nucleon self-energy, which are included in the SCGF approach yield a repulsive 
contribution the energy per nucleon, which increases with increasing nuclear density (dashed-dotted green line in 
Fig. [lj . This shifts the saturation point to a lower density and binding energy per nucleon so that also the saturation 
point obtained with SCGF is within the Coester band. 

The Hartree-Fock calculation using Vi ow k does not lead to a minimum in the energy versus density plot (see solid 
red line in Fig. [1]), as we have already mentioned before [2cj. |32T]. 

In order to reproduce the empirical saturation point of symmetric nuclear we have added an isoscalar interaction 
term as defined in eq. (f3"0")) choosing a value for a =0.5 and fitting the parameters to and t^. The results for these 
fitting parameters are listed in table |TT] and the corresponding energy versus density curves are displayed in Fig. [1] 
using the same line shape and color with and without adding the contact term but add an additional symbol to the 
lines displaying the results with inclusion of the contact term. For all three cases the fit yields an attractive two-body 
contact interaction and a repulsive £3 term. 

The results for the calculated saturation points in table Q] are supplemented by the corresponding values for the 
nuclear compressibility modulus 



dp 2 



(31) 

P=Po 



This nuclear compressibility, which is calculated at the saturation density p , together with the increase of energy at 
large density displayed in Fig. [1] characterize the stiffness of the EoS of symmetric nuclear matter. Comparing the 
different approaches we find that the relativistic features included in the DDRMF approach lead to the stiffer EoS 
around the saturation density as well as at higher densities. The SCGF and the Vi ow k calculations yield rather similar 
results after the contact terms are included, which are a little bit softer than the DDRMF results and characterized 
by a compression modulus of 270 MeV and 258 MeV for SCGF and Vi ow k, respectively. At higher densities the results 
are also very close to those obtained for the Skyrme Hartree-Fock using SLy4. Note, however, that SLy4 yields a 
rather low value for K as compared to the SCGF and Viowk calculations. The softest EoS for symmetric matter 
among those approaches which fit the empirical saturation point is provided by the BHF approximation. 

It is obvious that these results for the EoS for BHF, SCGF and Vi ow k are rather sensible to the choice of the 
exponent a in the £3 term of the contact interaction. A larger value of a would lead to stiffer EoS. We have picked 
the value a = 0.5 to obtain results for the EoS, which are similar to those resulting from the empirical approaches. 
Note, that the choice of a does not have any effect on the results referring to proton neutron asymmetries, which is 
the main focus of this study. 

Table [Til also displays results for the symmetry energy 



(32) 
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FIG. 2: (Color online) Comparison of the symmetry energy as(p) as a function of density p as obtained from 
Skyrme SLy4, DDRMF, BHF, SCGF, and Vi owk approaches. 



evaluated for each approach at the corresponding saturation density po- The two phenomenological approaches SLy4 
and DDRMF yield results which are in the range of the experimental value of 32 ± 1 MeV. Also the BHF and SCGF 
approach lead to results which are rather close to the empirical value, if the contact term has been added. The 
BHF and SCGF calculations without the contact term lead to non-realistic values for ag(p ) since these values are 
calculated at the corresponding saturation densities, which are larger than the empirical saturation density. 

The symmetry energy calculated in the SCGF approach is slightly smaller than the one obtained from the BHF 
approximation. This is valid for all densities under consideration (see Fig.[2|). This difference can easily be explained: 
As we already mentioned above, the contribution of the hole-hole terms is repulsive, which leads to larger energies for 
SCGF as compared to BHF for all densities in symmetric nuclear matter (Fig. [T|) as well as in pure neutron matter 
(see Fig. [3]). Since, however, the contribution of ladder diagrams is larger in the proton-neutron interaction (due to 
the strong tensor terms in the 3 <Si — 3 D\ partial wave) than in the neutron-neutron interaction, this repulsive effect 
is stronger in symmetric nuclear matter than in neutron enriched matter. Therefore the symmetry energy calculated 
in SCGF is slightly smaller if the hole-hole terms are included in SCGF(50j. 

The symmetry energy at saturation density obtained from Vi ow k plus contact term is only about two third of the 
experimental value (see tabled and the value is significantly below the other two microscopic approaches also at higher 
densities. This can be understood from the following considerations: The Hartree-Fock calculations using Vi ow k do 
not account for the attractive contributions due to the NN ladder terms involving NN states with relative momenta 
below the cut-off A. This missing attraction is compensated by the fit of the contact interaction to the empirical 
saturation point of symmetric matter. While the contact interaction is chosen to be identical for proton-neutron and 
neutron-neutron interaction, the ladder terms are more attractive for the isospin T = partial waves (see above), i.e. 
the proton-neutron interaction. This leads to a significant underestimate for the symmetry energy at all densities. 
Note, however, that this failure of Viowk should disappear if Vi ow k would be employed in an appropriate many-body 
calculation beyond the mean field approximation. 

The symmetry energy rises as a function of density for all approaches considered. Note, however, that the two 
phenomenological approaches Skyrme Hartree-Fock using SLy4 and DDRMF provide rather different predictions 
at high densities although the symmetry energy at normal density is identical. The relativistic approach predicts 
symmetry energies for high densities, which are well above all those derived from the microscopic calculations, while 
the Skyrme interaction yields a symmetry energy which is even below the Viowk estimate at densities above four times 
saturation density. 

Rather similar features also observed, when we inspect the properties of nuclear matter in /3-equilibrium, neutralizing 
the charge of the protons by electrons, displayed in Fig. [5] The upper panel of this figure displays the proton abundance 
Y p — Z/A, which are to some extent related to the symmetry energy: large symmetry energy should correspond to 
large proton abundances. So the largest proton abundances are predicted within the DDRMF approach. Already at 
a density around 0.4 fm -3 Y p exceeds the about 10%, which implies that the direct URCA process could be enabled, 
which should be reflected in a fast cooling of a neutron star. 
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FIG. 3: (Color online) Energy per nucleon of pure neutron matter as a function of density as obtained from Skyrme 

SLy4, DDRMF, BHF, SCGF, and Vi owk approaches. 




FIG. 4: (Color online) Results for a system of infinite matter consisting of protons, neutrons and electrons in 
/3-equilibrium. The upper panel show the proton abundances and the lower panel displays the energy per nucleon as 
a function of density using the various approximation schemes discussed in the text. 
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FIG. 5: (Color online) Effective masses for protons (lines with symbols) and neutrons (lines without symbols) as 
obtained for nuclear matter in /3-equilibrium using Skyrme HF (SLy4), BHF, and Vi ow k approaches. 



The other extreme case is the prediction derived from SLy4. In this approach the proton abundance does not exceed 
a value of 6%. 

The Viowk and SCGF approaches lead to similar proton abundances at large densities. This demonstrates that the 
evaluation of the proton abundance in /3-equilibrium cannot directly be deduced from the symmetry energy, since 
the former observable is derived from proton and neutron energies at large asymmetries (Z << N), whereas the 
symmetry energy is calculated from the second derivative at N = Z (see eq. (|32|) . The BHF approach shows slightly 
lower values for Y p at high density, but the results are still in the same range as SCGF and V[ ow k- 

At low densities the Skyrme HF approach yields large proton fractions as compared to the results of the other 
calculations. Large proton fractions at low densities tend to enhance density inhomogeneities and thus favor the 
existence of a large variety of pasta structures. Therefore the Skyrme HF (Sly4) and the DDRMF approach, which 
have been explored in detail in 0, [lol ] , should favor the formation of pasta structures as compared to the microscopic 
approaches. 

Comparing the energies of matter in /3-equilibrium derived from the various approaches as a function of density 
(Fig. flower panel) we find the same trends as in the case of pure neutron matter displayed in Fig. [3l The absolute 
values are lower in the case /3-equilibrium (about 75% of the energies for neutron matter) . Furthermore the relative 
differences between the various approaches are smaller. While in the case of neutron matter the energy differences 
between the various predictions are as large as 50% of the mean value, the corresponding number for matter in 
/3-equilibrium is only around 25%. The approximation schemes leading to large energies for neutron matter also 
show large symmetry energies, which result into relatively large proton abundances and smaller energies for matter 
in /3-equilibrium. 

The equation of state of nuclear matter in /3-equilibrium is the main input to predict mass and radii of neutron 
stars. A stiffer equation of state supports a larger maximum mass and a lower central density. In addition a thicker 
crust is found for the stiffer equation of state[51|. 

Another important information for the evaluation of dynamical features of matter in neutron stars is the density 
of states, which can be characterized by an effective mass. The term effective mass is used in various connections in 
many-body physics. This includes the effective Dirac mass of the relativistic mean field approach m* D (see eq. (|15[) ), 
as well as effective masses, which express the non-locality of the self-energy in space and time, which corresponds to 
a momentum and energy dependence. The density of states, however, is related to the single-particle spectrum close 
to the Fermi energy, which in the case of nuclear matter can be parameterized in terms of an effective mass by the 
expression 



e{k) 



k 2 



2m* 



U . 



(33) 



Such effective masses for protons and neutrons determined for nuclear matter in /3-equilibrium are displayed in Fig. [5] 
as a function of density considering non-relativistic approximation schemes. 
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It is a general feature of all approaches considered that the effective masses for protons as well as neutrons decrease 
with increasing density. However, there is a striking difference between the phenomenological Skyrme approximation 
and the BHF and Vi ow k approach, which are based on realistic NN interactions: The effective mass for protons is 
smaller than the corresponding one for neutrons in neutron rich matter for the calculations using realistic interactions, 
while it is opposite applying the Skyrme parameterization. In fact, if we define the effective masses for protons m* 



and neutrons 
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(34) 



it turns out most of the Skyrme parameterizations yield an effective isovector mass m v , which is even larger than 
the bare nucleon mass M[5|], which implies that it is larger than the effective isoscalar mass m* s . This means that 
the effective mass for neutrons is smaller than the corresponding one for the protons in neutron rich matter (/ > 0). 
These Skyrme parameterizations leading to a large effective isovector mass are usually favored as they correspond 
within the mean- field approach to an enhancement factor k of the Thomas- Reiche-Kuhn sum-rule [il |52|. Attempts 
have been made to distinguish between effective masses, which describe the energies around the Fermi energy, and 
those characterizing the bulk spectrum by introducing a term, which leads to a surface peaking of the effective mass 
term in finite nuclei [531] - 

Non-relativistic descriptions of nuclear matter, which are based on realistic interactions yield an effective isovector 
mass m v which is smaller than the corresponding effective isoscalar mass, which leads to a larger effective mass 
for neutrons than for protons in neutron-rich matter (see Fig. [5]). In order to analyze this finding we inspect the 
dependence of the nucleon self-energy in the BHF approximation Y J f HF 1 defined in eq.([7]), as a function of energy u> 
and momentum fc of the nucleon considered. Following the discussion of Mahaux and Sartor [54j one can define the 
effective fc-mass 



and the effective E'-mass 



m k (k) 
M 



m E (w) 
M 



MdT,(k,u) 
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dk 
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The effective mass can then be calculated from the effective fc-mass and the effective E-mass by 

m*(fc) TOfc(fc) m* E (u> = e(fc)) 



M 



M 



AI 



(35) 



(36) 



(37) 



Results for the effective fc-mass and E — mass as obtained from BHF calculations for asymmetric nuclear matter at 
a density p = 0.17 fm -3 and a proton abundance of 25 % (I = 0.5) are displayed in Fig. [6] We notice that the 
effective fc-mass for the protons is significantly below the corresponding value for the neutrons at all momenta. Since 
the fc-masses tend to increase as a function of the nucleon momentum fc, the difference in the Fermi momenta for 



protons and neutrons enhance the difference m£ n {k En ) 



:,p( k Fp)- 



The effective fc-mass describes the non-locality of the BHF self-energy. This non-locality and thereby also these 
features of the effective fc-mass are rather independent on the realistic interaction used. Furthermore it turns out 
that the values for the fc-mass are essentially identical if one derives them from the nucleon BHF self-energy using 
the G- matrix or from the bare interaction V or from Vi ow k [4l| . This non- locality of the self-energy is dominated by 
Fock-exchange contribution originating from 7r-exchange. In neutron-rich matter this contribution leads to a stronger 
depletion for the proton mass than for the neutron mass[20l. [55|. 

The effective .E-mass, representing the non-locality of the self-energy in time, yield values larger than M for momenta 
around kp- Also in this case the deviation of to* from M is larger for protons than for neutrons. The effective iS-mass 
originates from the energy-dependence of the G-matrix and is due to the admixture of 2-particle 1-hole configurations 
to the single-particle states. One may also say that the effective i?-mass is due to correlations beyond mean field. It 
accounts to the coupling of vibrational modes. 

Anyway, the enhancement of the effective mass m*, which is due to the effective E-mass in ea.p7p is not strong 
enough to compensate the effects of the fc-mass. Therefore the final effective mass is below the bare mass M and the 
effective mass for neutrons remains larger than the corresponding one for protons. 
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FIG. 6: (Color online) Effective fc-mass ml(k) (solid lines) and effective i?-mass m* E (k) (dashed lines) for neutrons 

and protons (lines with symbol) as obtained from the BHF calculations for asymmetric nuclear matter at the 
density p — 0.17 fm -3 and a proton abundance of 25 %. The Fermi momenta for protons and neutrons are indicated 

by vertical dotted lines 



The effects of the Z£-mass are weaker for Vi ow k than for the G-matrix. This is obvious since Vi ow k only accounts for 
ladder diagrams included with particle-particle states above the cut-off, whereas the G-matrix includes all particle- 
particle states above the Fermi-momenta. This explains the lower effective masses obtained for Vi ow k than for the 
BHF approximation (see Fig. [5|). Results on effective masses obtained from SCGF are rather similar to the BHF 
results, therefore we do not discuss them here explicitly. 

We want to add that the results for the effective i5-mass can be rather different in finite nuclei than in nuclear 
matter. As it has already been mentioned above, this i?-mass is due to the admixture of vibrational modes, which are 
quite different in finite nuclei as compared to nuclear matter. This may explain the differences in the predictions for 
the effective isovector mass resulting from Skyrme interactions, which are based on fits of properties for finite nuclei, 
and those originating from realistic interactions. 

Finally we are going to address the results for effective masses as they originate from relativistic mean-field approx- 
imations. Here we have to distinguish between the values for the Dirac mass and the effective mass parameterizing 
the single-particle spectrum according to eq. (|33p . The Dirac mass m* D is defined in terms of the scalar part of the 
self-energy (see eq. (|I5[l ). In the case of the relativistic mean-field approximation the self-energies do not depend on 
energy or momentum and the scalar self-energies for protons and neutrons are solely due to the direct contributions of 
the scalar mesons a and 8. The Dirac masses obtained from DDRMF in nuclear matter at /3-equilibrium are displayed 
in Fig. [7J We see that the Dirac masses decrease with density but show larger values for the protons than for the 
neutrons. 

In order to compare these results for the Dirac mass with the corresponding non-relativistic effective mass, we have 
to compare the expressions for the single-particle energies ea. ()17|i and eq. (|33p and adjust the parameters in such a 
way that the expressions yield identical results and slopes as a function of k at the Fermi momentum. This leads to 

< = ^Jk 2 Fi + (m* m )\ (38) 

where the label i refers to the case of proton and neutron. Results for these non-relativistic masses are displayed 
in Fig. [7J by solid lines. We find that the enhancement of the Dirac mass by the corresponding Fermi momenta in 
eg. (|3"8"]) is significant for the neutrons in particular. This leads to the effect that at high densities the non-relativistic 
effective mass for neutrons gets larger than the corresponding mass for protons, a behavior which is opposite to the 
one observed for the Dirac masses. 

It is worth noting that this feature, the difference of the Dirac masses m* Dn —m* Dp is negative while the difference of 
the corresponding non-relativistic masses is positive has been observed before within the framework of Dirac Brucckner 
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FIG. 7: (Color online) Effective masses originating from the DDRMF calculations of nuclear matter in /3-equilibrium. 
The Dirac masses for protons and neutrons arc represented by the dashed lines, while solid lines are used to identify 
the effective masses according to eq. . The results for protons are shown in terms of lines with symbols. 

Hartree Fock (DBHF) even at small densities [35l l39l l56l. [57| . The parameterization of the DDRMF approach has been 
made to reproduce the bulk properties of the DBHF of [3JJ. However, this adjustment cannot account for details like 
the non- locality of the self-energies in DBHF. Therefore it does not reproduce such details as the effective isovector 
mass with good accuracy. 

IV. CONCLUSION 

Various approaches to the nuclear many-body problem have been investigated to explore their predictions for nuclear 
matter at high density and large proton - neutron asymmetries. Two of these approaches, the Skyrme Hartree-Fock 
and the Density Dependent Relativistic Mean Field approach are predominantly of phenomenological origin. Their 
parameters have been adjusted to reproduce data of finite nuclei. However, the parameters have been selected in such 
a way that also bulk properties of asymmetric nuclear matter derived from microscopic calculations are reproduced. 
The other three approaches are based on realistic NN interactions, which fit the NN scattering phase shifts. In these 
approximation schemes (Brueckner Hartree Fock BHF, Self-consistent Greens Function SCGF and Hartree Fock using 
a renormalized interaction Vi ow k) a isoscalar contact interaction has been added to reproduce the empirical saturation 
point of symmetric nuclear matter. 

These various approximation schemes lead to rather similar predictions for the energy per nucleon of symmetric and 
asymmetric nuclear matter at high densities. In detail one finds that the relativistic DDRMF leads to a rather stiff 
Equation of State (EoS) for symmetric matter while the BHF approach leads to a relatively soft EoS, a feature which 
is compensated within the microscopic framework by the repulsive features of the hole-hole ladders included in SCGF. 
The phenomenological approximation schemes DDRMF (Skyrme Hartree Fock) over (under) estimate the symmetry 
energy at high densities as compared to the microscopic approaches. The lack of long range (low energy) correlation 
effects in Vi ow k leads to a symmetry energy which is too small already at normal density. These features are also 
reflected in the study of nuclear matter in the /3-equilibrium and lead to moderate differences in the predictions for 
proton abundances and EoS. 

More significant differences are observed when we inspect details like the effective masses, in particular the isovector 
effective mass. In neutron-rich matter the microscopic approaches predict a positive difference between neutron and 
proton effective masses. This feature can be related to the non-locality of the self-energy induced by one-pion exchange 
term and is expressed in terms of an effective A:-mass. This feature may partly be compensated by the effects of 
vibrational excitation modes on the nucleon mean fields. The effects of such low-energy excitations might lead to 
different results in nuclear matter and finite nuclei. This could be study e.g. in many-body calculations employing 
Viowk , which account for the effects of vibrational modes explicitly. 

We also discuss the differences between effective Dirac masses and corresponding non-relativistic masses in neutron- 
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rich matter. While the difference between the neutron and proton Dirac masses is negative, the differences of the 
corresponding non-relativistic masses tend to get positive. Further studies of the non-localities in space and time 
of various components in the Dirac self-energy would be useful to explore the connections to the non-relativistic 
microscopic approaches more in detail. 
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